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Abstract 

We discuss the issue of how precisely we can measure the primordial power spec- 
trum by using future observations of 21 cm fluctuations and cosmic microwave back- 
ground (CMB). For this purpose, we investigate projected constraints on the quan- 
tities characterizing primordial power spectrum: the spectral index n^, its running 
as and even its higher order running (ig. We show that future 21 cm observations 
in combinations with CMB would accurately measure above mentioned observables 
of primordial power spectrum. We also discuss its implications to some explicit 
inflationary models. 



1 Introduction 



It is now widely believed that the Universe experienced an extremely rapid expanding 
phase, the inflation, in the early stage of its history, which is assumed to be driven by 
the vacuum energy of a scalar field, called infiaton. In addition that the inflaton gives an 
inflationary expanding era, it can also provide the source of primordial density fluctuations, 
whose nature can be utilized to differentiate inflationary models. Specifically, the spectral 
index n^, describing the scale dependence of the primordial power spectrum, and the 
tensor-to-scalar ratio r, characterizing the amplitude of the gravity waves, are usually 
discussed, which can be measured/constrained by cosmological observations such as cosmic 
microwave background (CMB), large scale structure and so on. In particular, in the near 
future, several experiments are expected to measure/constrain the tensor-to-scalar ratio 
with high sensitivity, which may enable us to grasp invaluable hints to understand the 
inflationary model realized in the nature. 

However, on the other hand, there exist inflation models which predict too small tensor- 
to-scalar ratio to be detected, at least, in near future experiments. Furthermore, primordial 
density fluctuations are not necessarily generated from the inflaton, but they can be pro- 
duced by a light scalar field other than the inflaton, such as in the curvaton model [IHS], 
modulated reheating scenario [IjlS] and so on. In particular, models of these kind have been 
attracting attention since they can predict large values of /nl, a measure of bispectrum 



which characterizes non-Gaussian nature of density fluctuation^^. In these models, the 
tensor-to-scalar ratio is generally very small, and thus one may not be able to detect any 
gravity wave signal if such a model is the correct one for generating density fluctuations. 
In such a case, we may need to pursue another quantity which is not related to the gravity 
waves. As such a quantity, the spectral index is widely discussed. However, there are 
many inflationary models which can predict almost the same value of n^, thus we need 
something more. 

In light of these considerations, it would be desirable to have more observables, in 
particular, which are irrelevant to the tensor mode. Needless to say, such a quantity would 
also be helpful even if we can have information from the gravity waves as well, and give 
complementary information. As such an observable, we consider the running of the spectral 
index, as and even a higher order running or "the running of the running" [5s- The former 
quantity, Og, has already been explored by many authors. However, the latter one, /3s, has 



not been investigated much in the literaturd^n. on which we especially focus in this paper 



In general, higher order runnings are expected to be very small compared to the spectral 
index ?7.<j — 1, thus it seems very difficult to actually measure such a quantity. However, 
as we will show in this paper, future experiments of 21 cm fluctuation can give precise 
measurements of as and /3s, which would be very useful to differentiate inflationary models. 



Although the standard inflation models generally predict /nl ^{^)^ current observations of CMB 
and large scale structure give the constraints on /nl for the local type as —3 < /nl < 77 (95% C.L.) from 
WMAP9 [SI and 25 < /nl < 117 (95% C.L.) from NRAO VLA Sky Survey [5]. 

See Refs. [Hllin] for examples of such works. 
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In fact, the issue of probing the primordial power spectrum with 21 cm experiments was 
discussed in Refs. [II1[12]. In this paper, we further extend the discussion including the 
running of the running (5s and study its expected constraints from future observations of 
21 cm fluctuations in combination with CMB. In addition to observational constraints, we 
also study the runnings ag and /3s as well as Us in some explicit models and discuss its 
testability and in what cases the information from and Ps will be useful to differentiate 
inflationary models. 

This paper is organized into four sections as follows. In the next section, we give 
a formalism to discuss the power spectrum of primordial density fluctuations, paying 
particular attention to its scale dependence such as ns,as and (3s. We also give the 
predictions for these variables in some explicit models and discuss in what models/cases 
the information of the runnings are useful. In Section |3], we describe our method to derive 
expected constraints from future observations of 21 cm fluctuations and CMB and present 
out results. The final section is devoted to summary of this paper. 

2 Inflationary parameters: Formalism 

As mentioned in the introduction, the purpose of this paper is to investigate what kind of 
information can be obtained by considering the scale dependence of primordial curvature 
perturbations including "higher order" ones. Before going to study expected constraints, 
here in this section, we give a formalism and formulas to discuss primordial curvature 
perturbations and its scale dependence. We also study some explicit inflationary models. 
Since primordial curvature perturbation can also be generated from a light scalar field 
other than the inflaton such as in the curvaton, modulated reheating scenarios and so on, 
we give the formulas for such cases as well. But first we start with a case of the standard 
inflation model. 

2.1 Standard inflation case 

For the standard inflation case, the primordial curvature power spectrum is given, in terms 
of the 6N formalism, by 

where Vs^f, = {H^/2tiY (1 + 2(3C — l)e — 2Cri) is the power spectrum for fluctuations of 
the inflaton field [T3l[T3]. For the standard inflation case, dN/d((>^: = —H^/((>^ with a dot 
representing derivative with respect to time. The subscript "*" indicates that the quantity 
is evaluated at the time of horizon exit k = at^Ht, and C = 2 — In 2 — 6 ^ 0.73 with h being 
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the Euler-Mascheroni constant, e and rj are slow-roll parameters which are defined as 



^ = -Mi, 

2 ^^\V 



(2) 



(3) 



where V is the potential for the infiaton and a subscript represents a derivative with 
respect to the infiaton field. For later use, here we also define higher order slow-roll 
parameters: 



Pl Y3 



(4) 



(5) 



where we have put the superscripts (2) and (3) to remind that and a^^^ are second 
and third order slow-roll quantities, respectively. 

We expand the power spectrum in terms of the logarithm of the wave number up to 
the 3rd order as 



(n, - 1) In 



1 , 2f k 
-as m I - — 

2 V fcrcf 



where we have defined 



^ref 



^rcf 

ns-l+ios ln(fc/A,Vef)+|/3s In^Cfc/fcrof) 



3! V«^rcf 



/3s 
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(6) 



(7) 



as is usually called "running" of Ug. Thus in this sense, (3s may be called as "the running 
of the running." 

By using the slow- roll parameters, Us, as and Ps are written as 



10 \ 2 

- 1 = -6e + 2?7 + ( -y + 24C j + -r]^ 



a. 



-24e^ + 16er] - 2^ 



(2) 



(2 + 16C)er/+0 + 2C^e^2), (9) 

(10) 



(3s 



-1926^ + 192e'^r] - ?,2erf + (-24e + 2r/)^(2) ^ 2a^^\ 



'ID 
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2.2 Case with a light scalar field 

When the primordial curvature perturbation originates to fluctuations of a light scalar field 
cr other than the inflaton, such as in the case of the curvaton and modulated reheating 
models, the expressions for the spectral index and its runnings are different from those 
in the standard inflation case. Here we assume that the energy density of such a light 
scalar field is subdominant during inflation, and hence it does not affect the inflationary 
dynamics. To make our discussion general, we consider the case where the inflaton and 



a light scalar field can be both responsible for density perturbation^^. Cases with the 
standard inflation or a light scalar field being purely responsible for density perturbations 
are given as limiting cases. 

Denoting a light field as a and assuming that the inflaton and the a field are uncorre- 
cted, the total curvature perturbation can be written as 

V^{k)=vf{k)+Vf{k), (12) 

where V^'^ and V^^"^ respectively represent the contribution from the inflaton and a fields. 
The inflaton part V^^^ is actually the one given in Eq. ([T]). For the contribution from a 
light scalar field, V^'' is given by 

Vt\k)=(^^yvs.{k), (13) 

where Vsa = {H^/2'jtY (1 — 2(1 — C)e — 2Crj„) is the power spectrum for fluctuations of a 
light scalar field a [H]. Here the slow- roll parameter 77^ for a field is defined similarly to 
those for the inflaton 

where we assume that the scalar potential is given by the sum of those for and a as 
V^cj)) + U{a). For later use, here we also define other higher order slow-roll parameters as 

^3^2)2 ' y^"") 

where a subscript a denotes the derivative with respect to a field. 



^■^ Models of this kind are called "mixed" model. Such models have been investigated in the context of 
the curvaton model [TSHT^ and modulated reheating scenario pO] . 
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The spectral index Ug and its runnings as and /3s in this case are given as follows: 



1 =S„ 



-6e + 2r] + 
-2e + 27]^ + 



10 

y 



+ 24C 1 + ^7^2 _ (2 + 16C)er/ + Q + 2C ) e^^) 
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(17) 



a, = H<^(-24e2 + 16er/ - 24'^) 
-4S0S,(4e + 2r/<^-2r/^)2, 



;i8) 



= H<^(-192e^ + 192e2r/ - 32er/2 + (-24e + 2r/)^(2) + 2a(^)) 

E^{-6Ae' + 56e'7/^ - Ser/J + 2Ae\ - Seri^rj^ - 12e4'^ - ^e^^J^ + 2//.^^'^ + 2af ), 
12H^H^(2e - 7/ + r/.)(8e2 - 6er/^ + 2er/^ + ^^'^ - ) 

+8S^S^(S<^-S,)(2e-r/ + r/,)^ (19) 

where we have defined the fraction of V^^^ and V^'^^ in the total power spectrum, respec- 
tively, as 



V. 



V, 



jyW _|_ p(o-) ■ 



(20) 



By definition, the sum of these quantities is unity, S, 



1. 



2.3 Some example 

Now we discuss some explicit models and their predictions for the spectral index Ug and 
the runnings ag and /3s. To provide an example of in what cases the information from 
the runnings is useful, here we consider models with small tensor-to-scalar ratio, in other 
words, the information from the (scalar-mode) power spectrum only available. We will 
argue below that, in some cases, "running of the running" can be very useful to discriminate 
infiationary models. For descriptive purpose, we consider three models and discuss their 
predictions on Ug, ag and f3g. 



Model I : Curvaton with quartic chaotic inflation model 

In this model, we assume that the curvaton field is totally responsible for the cos- 
mic density perturbations. However, even in the curvaton model, the infiationary 
expansion is driven by the infiaton field and we need to specify the infiation model. 
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or the inflaton potential, to compute the power spectrum. For the inflation model, 
we adopt the chaotic inflation model with a quartic potential: 



v{ct>) = A0^ 



(21) 



where A is a coupling parameter. In fact, the quartic chaotic inflation is already 
excluded by current cosmological observations [7] when the inflaton is completely 
responsible for density perturbations. However, when the curvaton generates pri- 
mordial density fluctuations, it is still viable, and even more worth mentioning, such 
an inflation model is preferable for this case [171119] since it can give the spectral 
index of Us ~ 0.96, which is close to the central value from observational constraints. 
We should also note here that, although the curvaton model has been attracting at- 
tention because it can generate large non-Gaussianity, but it is not necessarily large: 
the curvaton can also predict small values of /nl in some parameter range. Thus, 
non-Gaussianity may not give useful information in such a case even for this type of 
model. 

For the curvaton sector, we assume that its mass is very small so that we can set 
?7o- = in the calculation of Us and so on. Furthermore, as mentioned above, the 
curvaton is assumed to be totally responsible for the primordial power spectrum, 
and hence we also set So- = 1 (S^ = 0). To realize this situation, fluctuations from 
the inflaton should satisfy P^^*^"* ^ 2.4 x 10~^, which means that the inflation scale 
is low and the tensor-to-scalar ratio becomes very small. 

In Fig. [H we plot the values of rig, as and Ps in this model as well as those for Model 
II which will be discussed in the following. In the figure, the predictions are shown 
in the n^-a^ (left panel) and the ag-^s (right panel) planes. Here the number of 
e-folds at the horizon exit is varied as 40 < iVg < 60. 



Model II :Mutated hybrid inflation model 

For this model, we assume that the inflaton is fully responsible for density pertur- 
bations and the inflaton potential has the following form: 



which is called mutated hybrid inflation model [21] and q is assumed to be a positive 
integer. ^ is some energy scale which characterizes the model. 

In Fig. [H we show the predictions for n^, as and Ps in this model along with those 
for Model I. For illustration purposes, here we take g = 8 and = M^\. In the 
figure, A'^e is varied in the range iVg = 40 — 60 to plot the lines. We note that, in both 
models, the tensor-to-scalar ratio is small as r < 10~^. By looking at the prediction 




(22) 



6 



cm vaton (quartic chaotic) 
mutated hybrid 



-5.0x10' 



b -1.0x10" 



-1.5x10" 



-2.0x10" 



0.94 



0.96 



0.98 



" curvatoii (quartic chaotic) 
mutated hybrid 



-5.0x10' 



-1.0x10'- 



-1.5x10'- 



-2.0x10 



3.0x10'"* -2.0x10'"* -1.0x10'"* 



Figure 1: Predictions in the Us-as (left) and the as-^s (right) planes from curvaton model 
with chaotic inflation (^((^) = A0^) and mutated hybrid model "with g = 4 and /i = Mpi. 
The number of e-folds is varied as 40 < iVg < 60 in this figure. 
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Figure 2: Predictions in the n^-a^ (left) and the /S^-a^ (right) planes from Type III hilltop 
inflation models. From the top to bottom, r]Q is taken to be 0.04, 0.03, 0.02, and 0.01. 
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in the n^-as plane, one can see that these models give very similar predictions for 
Us and Ug. However, if we look at the as^Ps plane, their values seems separated 
compared to those in Ug-as plane, which shows an idea that a higher order running 
(is can give useful information to differentiate inflationary models. To state this 
more quantitatively, we need to know attainable constraints on these observables 
from future cosmological observations, which is the topic in the next section. 

Model III : Type-III hiilltop inflation models 

As a last example, we discuss a model in which both \as\ and can be large. In so- 
called Type-III hilltop inflation models |22l|23] (see also Refs. pH - [26] and references 
therein), we expect a large running and a large running of running towards small 
scales. The potential shape is represented by 

v(4')-v.[i + -^m{^yxj^,. (23) 

The field (j) rolls from the top of the hill towards the origin during the inflationary 
epoch. Intermediately, the inflation ends by a water-fall mechanism. In Fig. [21 we 
plot the predictions in the Ug-ag (left) and the /Ss-a^ (right) planes, in case of g = 6 
and Vq^^/M^i = 1.6 x 10~^. Because the tensor-to-scalar ratio is tiny with an order 
of r ~ 2 X 10~^° in this model, it would be difficult to detect it. However, with future 
precise measurements of as and 13 s, we will be able to discriminate this model from 
others. 
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3 Analysis 



Now we discuss expected observational constraints on the quantities such as the spectral 
index n<j, its runnings as and from future cosmological observations by using Fisher 
matrix formalism. For this purpose, we investigate a joint constraint from 21cm fluctua- 
tions and CMB. We first briefly describe the formalism in each below, then present our 
results on the future constraints. 



3.1 21cm power spectrum 

Here we briefly sketch and summarize how we can constrain primordial power spectrum 
by using 21 cm observations. For detailed descriptions, we refer the readers to, e.g., 
Refs. dUEMQ]. 

The 21 cm line from the hyperfine transition of neutral hydrogen can be observed as 
the differential brightness temperature relative to the CMB temperature Tcmb: 

where Aiq ^ 2.85 x 10~^^s~^ is the spontaneous decay rate of 21 cm transition and V21 
is the 21 cm line (rest) frequency. Ts{x) is the spin temperature, which is defined by 
ni/rio = 3 exp(— T21/T5') with uq and rii being respectively the number densities of singlet 
loS'1/2 and triplet 1151/2 states of neutral hydrogen atom. Here T21 = /ic//cbA2i is the 
temperature corresponding to 21 cm line with A21 being its wavelength. dv\\/dr is the 
gradient of the velocity along the line of sight. 

Now we are going to consider fluctuations of Ti,{x). By expanding the hydrogen number 
density hh and the ionization fraction Xj (we also sometimes use xh = 1 — a^j: the fraction 
of neutral hydrogen) as Unix) = nnil + S{x)) and Xi{x) = Xi{l + 6x{x)), we can rewrite 
Eq. ([24D as 

1 dv. 



Tt{x) = T, (1 - x,(l + 6,ix)) (1 + 5{x)) 1^1 - , (25) 

where we have assumed that ^ Tcmb since we focus on the epoch of reionization during 
which this condition is well satisfied. is the spatially averaged brightness temperature 
at redshift z and given by 

/Qbh^\ f 0.15 l + ^\'/' 
Tfe ~ 27 mK ^— ^ , 26 

and dvr/dr is the peculiar velocity along the line of sight. 

By denoting the fluctuation in Tb as 6Th{x) = Th{x) — xnTb, the 21 cm power spectrum 
P2i{k) in the /c-space is defined by 

{6T,*{k)Sn{k')) = {27rf 6%k - k')P2i{k). (27) 
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z 


Xh 


6^ 


Rxx 


(^xx 


""^xx 


"x5 


Rxs 


OixS 








[Mpc] 








[Mpc] 




9.2 


0.9 


0.208 


1.24 


-1.63 


0.38 


0.45 


0.56 


-0.4 


8.0 


0.7 


2.12 


1.63 


-0.1 


1.35 


1.47 


0.62 


0.46 


7.5 


0.5 


9.9 


1.3 


1.6 


2.3 


3.1 


0.58 


2.0 


7.0 


0.3 


77.0 


3.0 


4.5 


2.05 


8.2 


0.143 


28.0 



Table 1: Fiducial values for the parameters in Vxxik) and Vx&{k) (See Eqs. (132!) and 

m) m- 

By treating the peculiar velocity 5.^ = {dv^ / dr){l / aH) as a perturbation and using that 
its Fourier transform can be given by = —fi'^6{k) with fi = k ■ n being the cosine 

of the angle between the wave vector and the line of sight, the power spectrum can be 
written as 

P2i(fc) = P^o{k)+fi''P^2{k)+fi''P^4k), (28) 

where k = \k\ and 

P^o = Vs5 - 2Vx5 + Vxx, (29) 
P^2 = 2 {Vss - Vxs) , (30) 
Pf.^ = Vss- (31) 

Here Vss = T^x^PsSyV^s = f^XiXnPxS and V^x = T^x^P^x and Pss,PxS and P^x are the 
power spectra defined in the same manner as Eq. 0271) for S and S^- Since 6 represents 
fluctuations in the hydrogen number density, Pss traces that of matter which includes the 
information on primordial power spectrum. Pxs and Pxx can be neglected if we consider 
the era when the intergalactic medium (IGM) is completely neutral. However, after the 
reionization starts, in which we are interested, these two spectra contribute significantly. 
Although a rigorous evaluation of these power spectra may need some numerical simula- 
tions, here we adopt the treatment given in Ref. [11], where VxS and Vxx are assumed to 
have a specific form to match radiative transfer simulations of Refs. [3T|[32]. Their explicit 
forms are the followings: 

'Pxxik) = hi, [I + axx{kRxx) + [kRxxf] Vss{k), (32) 
Vxs{k) = hi, e---(^-^-)-('=^-)>5,(A:), (33) 

where hxx-, hxs, cixxi Ixx and Uxs are parameters which characterize the amplitudes and the 
shapes of the spectra. Rxx and R^s correspond to the effective size of the ionized bubbles. 
The values of these parameters we adopt in the analysis are listed in Table [TJ 

Now we are in the position to discuss the Fisher matrix for 21 cm observations. First of 
all, we note that experiments of 21 cm radiation do not directly measure the wave number 
k nor the power spectrum in /c-space P2i{k). Instead, an experiment measures the angular 
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location on the sky and the frequency which can be specified by the vector 



e = e^e, + OyCy + A/e,(= 0^ + A/e,). (34) 

Here the frequency is represented by its difference from the central redshift 2; of a given 
redshift bin. Then we can define the Fourier dual of as 

u = UxCx + UyCy + u\\ez (= it± + u\\e-^ . (35) 

Notice that, since u\\ is the Fourier dual of A/, it has the units of time. Assuming that 
the sky is fiat^J, we can linearize the relation r and 0. Denoting the vector perpendicular 
to the line of sight as r^, we have the relations 

®L = rJdA{z,), A/ = Ar||/y(^,) (36) 

where dA{z^) is the comoving angular diameter distance and y{z) = X2i{l + z)"^ / H (z) . Ary 
is the comoving distance intervals corresponding to the frequency intervals A/. Then, the 
relation between k and u can be written as 

u± = dAk±, U|| = ykii- (37) 

The power spectrum of 6Th in w-space can be defined in exactly the same manner as that 
for fc-space by replacing k with u in Eq. fl271) . By using the relation between k and u, the 
power spectra in each space are connected as 

= ^^MMi''^'*'*- '''' 

We use the w-space power spectrum in the following analysis. 

With the power spectrum P2i{u), the Fisher matrix is given by 

where 6P2i{u) is the error in the power spectrum measurements for a pixel u, and pi 
represents cosmological parameters. To be conservative, when we differentiate P21 (u) with 
respect to cosmological parameters, we fix Vss{k) in Eqs. ( l32l) and ( l33i) so that constraints 
only come from the Vss{k) terms in P^o, P^2, P^4. The error of the power spectrum 6P2i{u) 
comes from sample variance and experimental noise, and is given by 

5P2i{u) = -— . (40) 



'^'^ Even if we consider all-sky experiments, the flat-sky approximation can be valid as long as the data 
are analyzed in many small patches of the sky [llj . 



11 



Experiment N,: 



ant 



Ae{z = 8) L 



'mm 



max 



FoV 





[km] 



[deg^] 



SKA 
Omniscope 



7000 
10*^ 



10 
1 



TT 8.6^ 6.8-10 
2.06 X 10^ 6.8 - 10 



Table 2: Specifications for interferometers of 21 cm experiments adopted in tlie analysis. 
For Omniscope, we assume that the effective collecting area Af, and field of view are fixed. 

The first term on the right hand side gives the one from the sample variance. Nc = 
27rk±Ak±Ak\\V [z) / {27t)^ is the number of independent cells in an annulus summing over 
the azimuthal angle. Here V{z) = dA{z)'^y{z)B x FoV is the survey volume with B being 
the bandwidth and FoV= X'^/A^ is the field of view of the interferometer with A^. oc 
being the effective collecting area. The noise power spectrum, denoted as Pn{u±) in the 
above formula, is given by 



where T^ys = 280[(1 + z)/7.4:]'^-^ K is the system temperature [33] which is dominated 
by the sky temperature due to synchrotron radiation, to is the total observation time 
and n{u±) is the number density of the baseline, which depends on actual realization of 
antenna distributions of each experiment. 

In our analysis, we consider the redshift range z = 6.75 — 10.05, which we divide into 
4 bins: z = 6.75 - 7.25, 7.25 - 7.75, 7.75 - 8.25 and 8.25 - 10.05. For the wave number, 
we set fcminii = 27r/(i/i?) to avoid foreground contamination [29] and take kma.x = 2 Mpc~^ 
in order not to be affected by nonlinear effect which becomes important on k > /cmax- 

To obtain the future cosmological constraints from 21 cm experiments, we consider SKA 
[31] and Omniscope [32] whose specifications are shown in Table HJ In order to calculate 
number density of baseline n{uj_), we assume a realization of antenna distributions for 
both arrays in the same manner as Ref. [11]. For SKA, we take 16% of the total antennae 
distributed with a filled nucleus of radius 189 m, and these antennae are surrounded by 
4% of the total antennae, which has an antenna density profile of an inverse square of 
the radius. These antennae are surrounded by a further 30% of the total antennae in 
a uniform density annulus of outer radius 5 km. We ignore the measurements from the 
sparse distribution of the remaining 50% of the total antennae that are outside the annulus. 
For Omniscope, which is a future square-kilometer collecting area array optimized for 21 
cm tomography, we take all of antennae distributed with a filled nucleus. In addition, we 
assume an azimuthally symmetric distribution of the antenna in both arrays. 



As mentioned previously, although 21 cm experiments have strong power to probe the 
primordial power spectrum, especially, on small scales, observations of CMB greatly help 




(41) 



3.2 CMB 
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experiment 


frequency 


beam ^fwhm 


At 


Ap 




[GHz] 


[arcmin] 


[fiK arcmin] 


[/iK arcmin] 


Planck 


70 


14 


137 


195 




100 


9.5 


64.6 


104 




143 


7.1 


42.6 


80.9 


CMBpol 


70 


11 


2.96 


4.19 




100 


8 


2.29 


3.24 




150 


5 


2.21 


3.13 



Table 3: Specifications for Planck and CMBpol adopted in the analysis. For CMBpol, 
we assumed the mid-cost (EPIC-2m) mission and only used three frequency bands (70 
GHz, 100 GHz and 150 GHz) for a realistic foreground removal. In the both observations, 
/sky = 0.65 is assumed. 



to determine other cosmological parameters such as energy densities of dark matter, baryon 
and dark energy, and so on. To obtain the future constraints, we consider Planck [36] and 
CMBpol [37] whose specifications are summarized in Table |3l 
The Fisher information matrix for CMB is given by [30 ] [38 ] 139] 



F 



(CMB) 



(2/ 



/sky Trace 



dpi 



dp. 



(42) 



where C/ is a covariance matrix of CMB, and pi represents cosmological parameters. We 
take the maximum multipole to be /max = 3000 in the analysis. The covariance matrix C; 
is expressed as 

Ci = I Cr Cr + N[ I , (43) 







Cf 




where (X = TT, EE, TE) is an angular power spectrum of each unlensed CMB mode, 
Cf"^ is one of deflection angle of CMB due to weak gravitational lensing, C™ is cross 
correlation between temperature and deflection angle, A'^;^ {Y = T, P) is a noise power 
spectrum of temperature or polarization, and Nf is one of deflection angle. The noise 
power spectrum NX is given by. 



N!'{u) = Al.exp[l{l + l)aUi^)] 



(44) 



for a single frequency band, where Ay can be found in Table |3]and crb(z/) = 6'fwhm/ v81ii2 
characterizes the width of the beam. When multiple frequency bands are used, where 
is given by the sum of its inverse as 



1 



(45) 
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In addition, we estimate the noise power spectrum of deflection angle A^;"^ by assuming 
that we use lensing reconstruction with the quadratic estimator [ID], and compute it by 
using a public code FUTURCMB |41], which adopts the estimator. 

Equipped with these formalism, we can calculate the Fisher matrix to obtain projected 
constraints on the spectral index Hs and the runnings ag and /3s, which will be presented 
in the next subsection. 



3.3 Future constraints 

Now we present our results for projected constraints on cosmological parameters, pay- 
ing particular attention to ns,as and (3^ which characterize the scale-dependence of pri- 
mordial curvature perturbations. In addition to these parameters, we also vary the 
standard cosmological parameters. Thus we explorer an 8-dimensional parameter space: 
Qa, Qfyh"^, h, r, A^, ra^, a^; A, where and Q}, are energy densities of cosmological constant 
and baryon, respectively, h is the Hubble parameter in units of lOOkm/s/Mpc and r is 
the reionization optical depth. In the analysis, we assume a flat universe, and fix the 
Helium abundance to he Yp = 0.25. Neutrinos are treated as massless particles and its 
effective number is also fixed to be A^eff = 3.046. As,ns,as and Ps are parameters which 
characterize the primordial power spectrum as already mentioned, and we normalize its 
amplitude as As{kroi) = 'P((/cref)/(2.18 x 10~^) [^H3] . As a reference scale to parametrize 
these quantities, we take fcrcf = 0.05 Mpc~^ in the most analysis presented below. How- 
ever, we also discuss how this reference scale affects the determinations of Ug, and /3s at 
the final part of this section. Then we set the fiducial values of these parameters (except 
as, (3s) near the best fit of the WMAP 7-year results [H], and = 0, f3s = 0, so that 
(Qa, ^bh^, h, T, As, Us, as, (3s) = (0.729, 0.02227, 0.703, 0.085, 1, 0.966, 0, 0). 
The total Fisher matrix is given by 

^_^mcm)^^(^MB)^ (46) 

where i,j are subscript representing cosmological parameters. We report our results for 
several combinations of experiments. As already mentioned above, although 21 cm ex- 
periments enable us to probe the primordial power spectrum very precisely since they 
can measure fluctuations smaller scales compared to those of CMB, regarding the de- 
terminations of other cosmological parameters, CMB has more strong power. Since the 
scale-dependence of the power spectrum has degeneracies with other cosmological param- 
eters, reducing the uncertainties of such parameters is also important to precisely measure 
Us, as and /J^. 

In Fig. [31 we show projected constraints in the n^-a^ and a^-Z^^ planes where 2a 
limits are shown for the analysis from Planck, Planck -|-SKA, Planck+Omniscope and 
CMBpol, CMBpol+SKA, CMBpol-|-Omniscope. Uncertainties for each parameter ns,a 
and Ps are summarized in Table HI We note that the la uncertainties are reported in the 
table. Although CMB experiments can already give a very precise measurement of Us at 
better than 0{1%), when one includes the data from 21 cm fluctuations, especially with 



s 
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6ns 




5as 




Sf3s 




Planck 


4.57 X 10" 


-3 


7.84 X 10- 


-3 


1.13 X 10" 


-2 


Planck+ SKA 


2.97 X 10- 


-3 


4.06 X 10- 


-3 


4.67 X 10- 


-3 


Planck+ Omniscope 


5.73 X 10- 


-4 


1.06 X 10- 


-3 


7.38 X 10- 


-4 


CMBpol 


2.12 X 10- 


-3 


2.44 X 10- 


-3 


4.53 X 10- 


-3 


CMBpol+SKA 


1.87 X 10- 


-3 


2.33 X 10- 


-3 


3.59 X 10- 


-3 


CMBpol+Omniscope 


5.28 X 10- 


-4 


9.90 X 10- 


-4 


6.92 X 10- 


-4 



Table 4: la uncertainties for ns,as and f3s from various data sets. We take /cref = 
0.05 Mpc~^ to derive these constraints. 

Omniscope, the precision improves by even one order of magnitude as seen from the figure 
and the table. The same also holds for the runnings as and Ps- In particular, for the 
case with CMBpol+Omniscope, one can probe the runnings as and (3s with the precision 
down to O{10~^). As discussed in the previous section, some inflationary model predicts 
as = (9(10-^) and /3s = O{10~'^), thus such models can be tested with future experiments 
of CMB and 21 cm fluctuations by using the runnings of the power spectrum. 

Our results show that the scale-dependence of the primordial power spectrum can be 
well probed when one uses the data from 21 cm fluctuations, which is consistent with [TTl 
[T2] , and further demonstrate that the information of "higher order" scale-dependence such 
as the running ag and the running of the running Ps can also give significant information 
on models of primordial density fluctuations. 

Finally we discuss how the choice of /crcf affects the determination of the parameters 
Us, as and 13 s- In Fig- El expected constraints on the Ug-as and as-(3s planes are shown for 
several values of /Cj-ef for Planck+ Omniscope and CMBpol+Omniscope. la uncertainties 
for each parameter are summarized in Table O As seen from the figure, the choice of fcref 
affects the uncertainties for n^, and Ps by a factor of a few and also changes the direction 
of the degeneracy. From the viewpoint of determining the scale dependence parameters, 
the optimal reference scale would be the one which gives uncorrelated constraints among 
the parameters. However, by looking at Fig. [6l we can see that the reference scale giving 
uncorrelated measures for the parameter sets (n^, as) and (a^, (3s) are different. In addition, 
even just considering CMB or 21 cm experiment alone, the optimal scale seems to also 
depend on the specification of the experiments. But a general tendency is that the optimal 
scales is around from 0.05 Mpc~^ to 0.1 Mpc~^. Thus we have mainly used the reference 
scale krei = 0.05 Mpc~^ to show our results. 

4 Conclusion 

We have investigated how precisely one can measure the power spectrum of the curvature 
perturbation from future experiments of 21 cm fluctuations and CMB. In particular, we 
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^rcf 


0.002 


0.01 


0.05 


0.1 


0.2 


0.5 


[Mpc^^] 














6ns 


3.74 X 10-3 


2.60 X 10-3 


5.28 X 10-4 


4.06 X 10-4 


4.73 X 10-4 


3.35 X 10-4 


6as 


1.47 X 10"3 


1.89 X 10-3 


9.90 X 10-4 


5.42 X 10-4 


2.64 X 10-4 


6.89 X 10-4 


5(3s 


2.46 X 10"^ 


6.05 X 10-^ 


6.92 X 10-4 


6.94 X 10-4 


6.93 X 10-4 


6.85 X 10-4 



Table 5: Expected la uncertainties of n^, and /3s from CMBpol+Omniscope for several 
values of k^ei- 



have studied projected constraints on the parameters characterizing the scale-dependence 
of the power spectrum such as the spectral index Ug, its running and the running of 
the running Ps- Although the former two parameters have been well explored in various 
context in the literature, the latter one, Ps, has not been studied much in connection with 
cosmological probes. 

Although the gravity waves or the tensor mode can give significant information to the 
inflationary Universe once it is detected, there are many inflation models, such as small- 
field models, predicting too small tensor-to-scalar ratio. In addition, in models with a light 
scalar field such as the curvaton model, modulated reheating scenario and so on, which 
are of interest due to the possibilities of their giving large /nl, the tensor-to-scalar ratio 
also tends to be very small. If one of these models is realized in the nature, it would be 
very difficult to detect the signature from the gravity waves. However, even in that case, 
"higher order" scale dependence of (scalar) curvature perturbations would help to probe 
the inflationary model, which we show quantitatively in this paper. We have discussed 
some explicit models where higher order running Ps would be very useful to differentiate 
models. Needless to say, even when the tensor modes are detected, the runnings can give 
extra valuable information on models of primordial fluctuations. 

We have obtained expected constraints on such a higher order running /3j, as well 
as Us and by using observations of 21 cm fluctuation, in combination with CMB. 
Since the power spectrum of 21 cm fluctuations can probe cosmic density fluctuations on 
smaller scales than those observed in CMB, one can obtain severe constraints even for 
the running of running parameter Ps- In particular, when one considers the combina- 
tion of CMBpol+Omniscope, we can probe the running parameters with the precision of 
6as = O{10~^) and 6Ps = 0(10-4), which would give useful information to discriminate 
inflationary models. 

Although current cosmological observations are already so precise that some of infla- 
tionary models have been excluded, there are still many possibilities allowed and thus we 
need to go further to understand the early Universe more. In particular, to see the details 
and differentiate models well, it is preferable to have yet another observables other than 
commonly used one. As such a quantity, we considered a higher order running Ps in this 
paper, which can be well probed by future cosmological observations such as from 21 cm 
and CMB. In the view that cosmological data will be much more precise in the future, the 
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research along this hne would provide us a lot of insight on the inflationary Universe and 
even the origin of the Universe. 
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Figure 3: Projected constraints from Planck, Planck+SKA, Planck+Omniscope (top pan- 
els) and CMBpol, CMBpol+SKA, CMBpol+Omniscope (bottom panels) in the n^-a^ 
(left) and ag-^s (right) planes. 



21 



0.1 



0.05 



a 



-0.05 



-0.1. 



k_ref= 0.002 MpcA-1 kjef = 0.1 MpcA-1 

.k_ref = 0.01 MpcA-1 k_ref = 0.2 MpcA-1 

■.k_r>=.0.05 MpcA-1 k_ref=0,5MpcA-l 




Planck 



0.85 



0.9 



0.95 



1.05 



1.1 



0.04 



0.02 



C3 



-0.02 



-0.04 



k_ref = 0.002 MpcA-l k_ref = 0.1 MpcA-1 

k_ref = 0.01 Mpc"-1 kjef = 0.2 MpcA-1 

.k_ref=0.05 Mpc-^-l k_ref = 5 MpcA-1 



CMBPol 



0.9 



ns 



0.95 1 

ns 



1.05 



0.1 



0.05 



a - 



-0.05 - 



-0.1 



k_ref = 0.002 MpcA-1 k_ref = 0,1 Mpc^-l 

k_ref=0.01 MpcA-1 k_ref = 02 MpcA-1 

k_ref=0.05 MpcA 1 kjef = 0.5 Mpc^-l 




Planck 



-0.1 



-0.05 





/3s 



0.05 



0.1 



0.04 



0.02 



a 



-0.02 



-0.04 



k_ref = 0.002 MpcA-l k_ref = 0.1 Mpc^-l 

k_ref = 0,01 Mpc^-l kjef = 0.2 MpcA-1 

k_ref = 0.05 MpcA-1 — kjef = 0.5 MpcA-1 




IMBPol 



-0.04 -0.02 





(is 



0.02 0.04 



Figure 4: Projected constraints in the n^-as (top panels) and ag-^s (bottom panels) planes 
for several values of k^.^ from Planck(left) and CMBpol(right). 
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Figure 5: Same as Fig. HI but from Planck+SKA(left) and CMBpol+SKA(right). 
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Figure 6: Same as Fig. HI but from Planck+Omniscope (left) and CMBpol+Omniscope 
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